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Abstract. We describe the gated bilayer graphene system when it is subjected 
to intense terahertz frequency electromagnetic radiation. We examine the electron 
band structure and density of states via exact diagonalization methods within 
Floquet theory. We find that dynamical states are induced which lead to 
modification of the band structure. We first examine the situation where there 
is no external magnetic field. In the unbiased case, dynamical gaps appear in 
the spectrum which manifest as dips in the density of states. For finite interlayer 
bias (where a static gap is present in the band structure of unirradiated bilayer 
graphene), dynamical states may be induced in the static gap. These states can 
show a high degree of valley polarization. When the system is placed in a strong 
magnetic field, the radiation induces coupling between the Landau levels which 
allows dynamical levels to exist. For strong fields, this means the Landau levels 
arc smeared to form a near-continuum of states. 
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1. Introduction 

Graphene [I] and its bilayer [2] have attracted much attention recently due to the 
novel fundamental physics that they display and huge scope that they have for 
device applications (3J [3] - In particular, electrons in monolayer graphene (whose 
low-energy quasiparticles mimic chiral Dirac fermions with linear dispersion) have 
displayed relativistic-like phenomena, including Klein tunneling [6] and the half- 
integer quantum Hall effect [7J 15] ■ Electrons in bilayer graphene show properties which 
are hybrid between the monolayer and traditional two-dimensional semiconductor 
systems, since the low-energy quasiparticles are chiral, but the inter-layer coupling 
induces an effective mass and corresponding quadratic energy dispersion. A fully- 
tunable gap can be opened at the charge neutrality point by application of a transverse 
electric field (either by gating [9] or by doping [2]), a feature which is unique to this 
system. 

Monolayer graphene which is irradiated by monochromatic, high-intensity laser 
light has been studied theoretically for the bulk [TU1 EJ HH E3 HH [T5] , in nanoribbons 
|16H17j . and in np-junctions |18II19) . Experimental investigations of few-layer graphene 
have also been carried out [20]. However, the transport properties of bilayer graphene 
which is irradiated by intense laser light have not been considered in much detail. 
Ryzhii et al. proposed that a phototransistor could be implemented using bilayer 
graphene |21j . Wright et al. have found a large peak conductance in the terahertz 
and far infra-red frequency range for bilayer graphene nanoribbons [22 ], and we have 
previously suggested that valley-polarized electrons can be produced in gapped bilayer 
graphene [23] . 

In this paper, we give a comprehensive description of the band structure and 
density of states of irradiated bilayer graphene, both in zero external magnetic field and 
in a strong field. This is of fundamental physical importance, but also has application 
in the realm of devices and technology because of the growing consensus that graphene 
and its bilayer have vast potential in the fields of optoelectronics and photonics [4], 
and in the design of new electronic devices such as ambipolar transistors [24]. Also, the 
spin- like degrees of freedom in graphene (such as the lattice pseudospin and the valley) 
may allow for electronic implementations of the ideas of spintronics |25] which have 
been discussed in the literature [H [26] [2JJ 023 US] . In all of these thorough 
understanding of the basic properties of irradiated bilayer graphene is an essential 
building block for design and application of devices. In particular, this present work 
focuses on the generation of valley-polarized states which may be used as a filter for 
the generation of valley-polarized currents. This is an essential step in the realization 
of valley-tronics devices. 

We briefly outline the structure of our paper. In Section [2] we describe the 
theoretical framework which we employ for both the zero-field and strong-field cases. 
Then we present and discuss the results of our calculations in Section [3] before 
summarizing and placing our work in the context of valley-tronic devices in Section 
[5] Various important formulae and derivations are collected in the Appendix. 

2. Theoretical framework 

The Hamiltonian of irradiated bilayer graphene is written as 
U{t) = H + Hu + h(t) 
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where H is the continuum limit of the tight binding Hamiltonian, Hjj describes the 
external electrostatic field due to a top gate or dopants, and h(t) is the time-dependent 
part which depends only on the irradiating field. We take the nearest-neighbour 
approximation of this Hq, but in principle any combination of hopping terms can be 
included by computing the eigenfunctions of the static Hamiltonian numerically. The 
static part of the Hamiltonian Hq + Hjj determines single particle wave functions (px , 
which span the spatial part of the Hilbert space of solutions ofH(t) and which we shall 
use as a basis for the time-dependent solutions ^(t) of H(t). The energy spectrum 
and wave functions of the static Hamiltonian are presented in |Appcndix A| both in 
zero external magnetic field and when a strong magnetic field quantizes the electron 
motion into Landau levels. The solutions of the time-dependent part are found via 
the Floquet theorem, which we describe below. 

The time-dependent term in the Hamiltonian is not known a priori, so we describe 
the process by which it is derived. The irradiating field is represented as a time- 
dependent vector potential which is introduced to the static Hamiltonian via Peierl's 
substitution p — > p + eA, where e > is the magnitude of the electron charge. We 
write the vector potential as A(t) = [A x (t), A y (t)}, and substituting into the static 
Hamiltonian in Equation (|A.1[) to find 

H^^( Mt) l Mt) ^-^) (!) 

where a x is the first Pauli matrix. The vector potential encodes all information about 
the radiation, such as the polarization of the light, strength of the field, and frequency 
of the radiation. In the case of linearly polarized light, this field can be represented 
as Aii n — Acosf^i [cos 0, sin 0] where 9 is the angle of polarization in the plane of 
the graphene and Qa is the frequency. The strength of the vector field is given by 
the parameter A = F/Q.A where F = |E| is strength of the electric field. However, 
when the light is circularly polarized, there is a phase difference between the x and y 
components of the field so that 

A±= A [cos Q A t, ± sin Q A t] (2) 

where the plus or minus sign refer to right-handed and left-handed polarizations, 
respectively. Full details of the time-dependent Hamiltonians in each case are given 
in |Appendix B| 

It is not possible to develop exact solutions to the eigenvalue equation for H(t). 
However, since the pertubation h(t) is periodic, we may employ the Floquet theorem 
[50] to write the full time-dependent wave function ^(r, t) as the product of a periodic 
function $(r, i) (with period to = 2-k/VLa) and a time evolution function. The 
'temporal Brillouin zone' is the interval < t < to, and the wave function is 

<&(r,i) = e- fe */ R $(r,i) (3) 

where + t ) = $(r,i). The scalar constant e is called the quasienergy. This 

theorem is the temporal analogue of the Bloch theorem, so that the quasienergy 
is equivalent to the quasimomentum, and the time period is equivalent to the 
lattice constant of the reciprocal lattice. The periodic wave functions $ are called 
'Floquet states'. Substituting Equation ([3j into the time-dependent Schrodinger 
equation yields %(t)^(t) — ih-^^(t) = e<fr(t), an eigenvalue equation for the operator 
F{t) = Hit) - ih§- t with F{t)^{t) = 0. 

In order to solve the Schrodinger equation for J-(t) and we consider an 

expanded Hilbert space 1Z ® T of square-integrable functions of space and functions 
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of time with period to (see Ref. [3T] for a full description). In this space, the scalar 
product is defined as the regular spatial scalar product with the average over one 
period: 



1 



to 



((n\m)) = — / (n\m) dt. 
to Jo 

where n and m label arbitrary states in 1Z <£> T . Since the wave functions are periodic 
and the eigenfunctions of the static Hamiltonian Hq + Hu form a complete set for the 
spatial coordinates, we can write the wave functions as 

S» TO (r,t) =J2c?x*jx(r,t) (4) 

j,x 

where m labels the Floquet state and the Hilbert space 1Z ® T is spanned by the 
infinite set of functions 

i J x(r,() = e'^VxW 

such that j G {...,—2,-1,0,1,2,...}, and X labels the eigenstates of Ho + Hjj. 
The label X contains all appropriate single-particle quantum numbers, but its exact 
composition depends on whether there is a magnetic field present in the system. The 
operator in Equation ([1} does not couple states with different momenta, different 
spin, or which are in different valleys. Therefore the Floquet states which result 
from diagonalization of the ^-dependent Hamiltonian retain these three quantities as 
good quantum numbers. The Hamiltonian can be written as a matrix by computing 
the matrix elements of T(t) over these states. This yields an infinite dimensional 
matrix which can be truncated for a sufficiently large number of terms in the Fourier 
expansion and numerically diagonalized to give quasienergies and wave functions to 
arbitrary precision. These matrix elements are discussed in |Appcndix B| 

We now introduce the two-time formalism which we use to compute the Green's 
function and hence the density of states (DoS) in irradiated bilayer graphene. In this 
formalism, the time associated with the expanded Hilbert space (which was previously 
labelled 't' but which we shall call £ from now on) is separated from the evolution 
of the system such that the full time-dependent solution of the Schrddinger equation 

¥(r,t) is ism 

*(r,t) = tf(r,C,t) _. (5) 
The two-time wave function is then defined to be 



¥(r,C,t)=expf-^F(C)(t-*0) *(r,C,0 



where T is the Floquet operator introduced earlier. Full time-dependent solutions are 
given by the limiting procedure in Equation ([5]), but we shall generally be interested 
in the dynamics of the system on timescales much longer than to so we shall instead 
take the time average with respect to the field. 

In the Matsubara formalism, we utilize the grand canonical ensemble, and define 
the associated energy scale k = e — fi where /i is the chemical potential. The operator 
for this energy is K. = J 7 — fj,N. The imaginary time r = it is defined and the evolution 
of the field operators associated with the Floquet states is given by (the coordinate 
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dependence is implicit) 



^(C,T)=5>-^<MC)a. 



(6) 
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(7) 



n 



where the index n runs over all Floquet states. The Matsubara Green's function is 



The operator T T is the r-ordering operator, f2 is the thermodynamic potential and 
serves as the normalizing factor for the thermodynamic average, and /3 = l/(/csT). 
The Fourier transform of this Green's function is the quantity from which the density 
of states can be calculated. The Fourier transform is 



which includes an averaging procedure over the period of the fast oscillation associated 
with the radiation. The retarded Green's function G R can then be found by carrying 
out the analytic continuation iu n — > ui + iS and the density of states can be extracted 
from this Green's function in the standard way: 



where tr denotes the summation of the diagonal elements of G R which is a 4 x 4 matrix 
in the sublattice space. 

3. Results 

3.1. Zero magnetic field 

In the case of zero external magnetic field, the single particle quantum states in the 
system are of plane wave nature and are characterized by a two-dimensional wave 
vector k. In addition, there are two (real) spins cr, two inequivalent valleys £, and four 
bands within each valley (labeled by the conduction or valence band v and the high- 
energy or split branch b) . The wave functions associated with these states are given in 
| Appendix A.l| The energy spectrum is found by substituting the matrix elements 
of the radiation operator into the Hamiltonian and numerically diagonalizing the 
resulting matrix. The spectrum is shown in Figured! Panel (a) shows the unirradiated 
low energy bands for comparison with the other three plots. In zero magnetic field, the 
effective coupling parameter is vpeF /ftQ 2 A which implies that the strongest coupling 
occurs for smaller frequencies. This is illustrated by panels (b) and (c) which show 
the effects of an irradiating field with Ha = 2 THz and 5 THz respectively with a 
field strength of F — 5kVcm _1 . The mixing between different Fourier harmonics is 
stronger for the lower frequency, but the dynamic gap opened by the mixing is larger 
for the higher frequency. The dynamic gaps only occur when states from opposite 
bands mix. States from the same band run parallel to each other (seperated in energy 
by MIa) and therefore cannot cross. There is no gap at k — because the distribution 



Q (rC, r'C, t-t') = — Tr e^^T^C, r)^ (C, r') 





(9) 



p 
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Figure 1. Energy spectrum of irradiated bilayer graphene in zero external 
magnetic field under right-handed circularly polarized light, (a) The unirradiated 
spectrum; (b) Q A = 2THz, F = SkVcm" 1 ; (c) U = 100 meV, Q A = 2THz, 
F = SkVcm" 1 ; (d) Q, A = 5THz, F = SkVcm" 1 . 
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Figure 2. Valley asymmetry when interlayer bias potential is present. In these 
plots F = 5kVcm-\ Q A = 2THz, and U = 30meV. 



of the wavefunction across the four lattice sites forbids coupling for small momenta 
(see |Appendix A[ ) . 

When a gap is introduced to the spectrum by doping one side of the graphene 
bilayer or by electrostatic gating, the effect of the radiation is markedly different. This 
is because the gap means that consecutive Fourier harmonics from opposite bands do 
not cross, and therefore the spectral weight is spread only between Floquet states 
which originate from the same band. Therefore no anticrossings appear, although the 
effect of a broadening of the band will be seen. This is illustrated in panel (d) of Figure 
[1] However, as Figure [5] shows, the coupling to right-handed circularly polarized light 
is stronger in the K valley than in the K' valley. This is due to the redistribution 
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of the wave function among the four lattice sites as a result of the interlayer bias 
potential which generates the gap. If the direction of the potential or the orientation 
of the irradiating field are changed, then the coupling becomes stronger in the K' 
valley instead. If linearly polarized light is used, the response of electrons in the two 
valleys are identical since linearly polarized radiation can be represented as the sum 
of two circularly polarized components. 

The Floquet wave functions which are calculated numerically can now be used 
in Equations ^ and © to compute the Green's function and hence the density of 
states. The Green's function is 



where n' labels the discrete set of Floquet states with wave vector p and x labels the 
unirradiated basis states with wave vector p. The density of states can be calculated by 
substituting this expression into Equation ([9]) and numerically evaluating the integral 
over momentum: 



where A is the momentum cut-off determined by requiring the summation over all 
states yields the correct electron density at half-filling. Plots of this function are 
shown in Figure [3l In (a), the interlayer bias potential is zero so that the graphene 
does not have a static gap. The density of states for unirradiated graphene is constant 
in this case, as predicted by straightforward analysis of the band structure. When 
the graphene is irradiated, dynamic gaps open at and are clearly visible for 

stronger fields. The responses of electrons in the two valleys are identical. In (b), 
the same system parameters are used, except that now a small static gap is present 
(U = lOmcV = U/2). In the unirradiated case (thin solid line), the static gap 
manifests as a region with near-zero DoS for small energies (fiw < 5meV). For 
hu > 5mcV, the dynamic gaps are still present. However, for hw < 5meV, a finite DoS 
is present under strong radiation. This is due to the dynamical states that are induced 
by the radiation. In (c), the gap is wider (U = 20meV) and now the conduction and 
valence bands are seperated to an extent where significant coupling between them is 
not permitted. However, states within the same band do still couple causing the band 
edge to be significantly smeared and for electron density to be present in the static 
gap. 

3.2. Quantizing magnetic field 

When a strong magnetic field is present in the system, the motion of the electrons is 
quantized into Landau levels. Using the Landau guage Ab = [0,Bx,0], the discrete 
single particle quantum numbers are the same as in the zero field case, but the 
momentum is continuous in the y direction and discrete in the x direction. This 
is due to the gauge field which breaks the translational symmetry. Therefore, the sum 
over the two-dimentional momenta may be split into two separate one-dimensional 
sums, one discrete (over k x ) and the other represented by an integral with periodic 
boundary conditions (over k y ). 

Figure |4] shows the evolution of the Landau level spectrum with increasing 
intensity of linearly polarized incident radiation. At F — 0, the standard bilayer 
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Figure 3. Density of states in bilayer graphene in zero magnetic field, (a) 
For ungappcd bilayer graphene, the valley degeneracy is intact and gaps open at 
intervals of £2^/2 ~ 4. lmeV. (b) U = lOmeV. A significant electron density is 
induced in the spectral gap by the radiation, (c) U = 20 meV. The spectral gap 
is too wide to allow significant coupling of electron and hole states. 
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Figure 4. Evolution of energy levels with the strength of a linearly polarized 
irradiating field. The line thickness indicates the strength of the mode, Q = 1 THz, 
U = 20 meV, and B = 5 T. The spectrum of levels is symmetric about zero energy 
in this case. 



Landau level spectrum is evident. For weak field, (F < lkVcm 1 ) the mixing of 
dynamical states is small and the original Landau level spectrum is recognizable. 
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Figure 5. Density of states for B = 12T, U A = 4THz, U = 20meV, and with 
F = 5kVcm~ 1 in each valley. The radiation shifts the energy of the Landau 
level, and coupling properties are different in the two bands. 



For strong coupling (F > SkVcm^ 1 ), the Landau level spectrum is replaced by a 
near-continuum of levels, each with rather small weight, and there are the beginings 
of states visible in the gapped region. Notice, however, that the two 5-states barely 
change their intensity showing that they are only very weakly coupled to the radiation 
due to their unique distribution of wave function weight between the four sublattice 
sites. 

The DoS in this case is evaluated using the same steps as in the zero magnetic 
field situation. The analytical expression which we must evaluate for the DoS is 

P = 77T - 77 72 77 \ c ox I 

7rAs 1C ~ Ka > + 6 x 

This function is plotted in Figure [5] for bilayer graphene with a small inter-layer bias 
U = 20meV for left-handed, right-handed and linear polarization of the incident light. 
We clearly see the radiation-induced dynamical bands as expected. The first feature 
of these plots which requires discussion is the difference between the behavior in the 
conduction and valence bands for the two circular polarizations. This occurs because 
coupling between the radiation and the light and the electrons depends critically on 
the wave function components. The factors like £5 ± e mean that the weight on each 
lattice component are different in the two bands, implying that the field will affect 
these electrons differently. When strong coupling occurs, the Landau level spectrum 
is almost washed out by the dynamical bands, but when coupling is weak, the states 
are still clearly discernable. There is also a slight shift in the energy of the Landau 
levels in each valley, which is caused by the broken valley degeneracy in the single- 
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Energy [meV] 

Figure 6. Valley polarization for zero magnetic field, a) Polarization for right- 
handed (solid line) and left-handed (dashed line) radiation with Q A = 2 THz, 
F = lkVcm -1 , and U = 20meV. b) Demonstration that polarization of almost 
unity can be achieved under specific conditions. Solid line has il A = 4 THz, 
F = SkVcm" 1 and U = lOmeV; Dashed line is Q A = 6THz, F = 8kVcm-\ 
and U = 20meV. c) Polarization features caused by the dynamic spectral gaps at 
u) = become more pronouned with increasing field stength. In this panel, 

n A = 6 THz and U = 20meV. 



particle spectrum [see Equation (|A.3j) ]. The linear polarization shows this as a uniform 
downward shift of the Landau levels in the K valley relative to the K' valley. 

4. Conclusions and summary 

We have shown that under certain conditions, the coupling of electrons in biased 
bilayer graphene to the external electromagnetic field is not the same in either valley. 
In order to demonstrate how this might be applicable in the design of a device, we 
discuss how the existence of parameter ranges where there is significant DoS in one 
valley but not in the other allows us to consider the possiblity of generating valley- 
polarized electrons [23) . If a current is incident on a region of irradiated bilayer 
graphene where the radiation and inter-layer bias are tuned such that states only 
exist in one valley, then incident current in the valley where there are no states will 
not be able to traverse the irradiated region. This region can then act as a filter for 
the electron valley, leading to the possibility of switching devices and 'valleytronic' 
applications [23] . We define the polarization to be 

„ _ PK ~ PK> 
PK + pK> 
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so that V = +1 implies that all electrons are in the K valley and V = — 1 means 
that all are in the K' valley. Figure [5] shows the polarization for several different 
scenarios. Firstly, in panel a) we demonstrate that the left-handed and right-handed 
orientations of light induce polarization in opposite valleys. In panel b) we demonstrate 
that polarization can approach unity for specific values of the radiation parameters 
and bias potential. This is important if highly-polarized electron currents are to be 
produced. Lastly, in panel c) we show that the polarization induced near the dynamical 
gaps becomes more pronounced and exists for a wider energy range when the field 
strength is increased. This is in contrast to polarization features caused by additional 
dynamical states appearing in the static gap, which exist only for small parameter 
ranges. The features associated with the dynamical gaps also become stronger as the 
frequency is increased. 

In the case where a strong magnetic field is present, valley polarization may also 
be generated, but only as a result of the shift in energy of the Landau levels created 
by the radiation. 

This work was supported by the Canada Research Chairs programme and the 
NSERC Discovery grant. 



Appendix A. Single particle description 

We take the basis {tpA,,'fpB u ,^A u ,^B l } with £ = +1 in the K valley and the basis 
{tpB^TtpAntpBnipA^} with £ = — 1 in the K' valley so that the single particle nearest 
neighbour tight binding Hamiltonian for unbiased bilayer graphene can be written as 




















7i 





7i 




/ 



(A.1) 



where n — p x + ip y is the linear expansion of the transfer integral in the tight binding 
formalism, 71 is the energy of the inter-layer dimer bond, and vf is the Fermi velocity. 
Similarly, the inter-layer bias potential which sets the energy of lattice sites in the 
upper layer at U/2 and sites in the lower layer at — U/2 causes a gap of magnitude U 
to open at the K points, and is described in [33J 

/ (U/2 \ 

-£U/2 

-(U/2 

\ (U/2 ) 



Hj, 



(A.2) 



Appendix A.l. Zero magnetic field 

In the zero field case, the 7r operators are constructed from the usual single-particle 
momentum operators it — —ihd x — hd y . The energy spectrum associated with the 
Hamiltonian H + Hjj is 




h 2 v 2 F k 2 + bJ^ + h?v F k*{U 2 + 7l 2 ) 
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The quantum numbers v = ±1 and b = ±1 label the band and branch respectively. 
Henceforth, we denote the band and valley indices by the label a — {£ Q ,^a,&a} to 
shorten the notation. The wave functions are given by 



<Pk,a = Ck,c 



1 



V 



-Vi ^ ^ 



k.<v 



/ 



/k,a 



h 2 v 2 F k 2 



h 2 v 2 F k 2 



(£ Q s + £ k>a ) 2 



(£ Q S + £ k , Q ) 2 
+ 1 



(£«£ — E^ a ) 2 



(£q£ + E'k.q)^ 

7? 



/k,a 



where H = U/2 and # k is the angle of the wave vector k in the graphene plane. 



Appendix A. 2. Quantizing magnetic field 

In this case, the 7r operators are constructed from the guage-invariant momentum in a 
magnetic field found by making the minimal coupling substitution in the momentum 
operator so that p — > p + eA. We label the Landau levels with a set of quantum 
numbers a = {n a , u a , £ a , q a } where n > is the Landau level index, v = +1(— 1) in the 
conduction (valence) band, £ is the valley, and q is the x component of the wave vector 
which defines the guiding centre coordinate in the Landau gauge and we assume that 
all Landau levels are in the low-energy bands. Then, the Landau level energies E a 
are found by solving a polynomial equation derived using the Landau level operators 



n Xn = -iH^-Xn-i and ^ X 



.hJ2(n+l) 
l^-^Z Xn+V 



E 2 a )=0 (A.3) 



Ab A.n-L » An — ' x B 

[2(n a + 1) - (Hf a - E a ) 2 ] [2n a - (S |o + E a ) 2 ] - 7l 2 (~ 2 a 

where the energies are measured in units of Hvf/^b and = ££7/2. The wave 
function associated with each Landau level is defined by the index of the Landau 
function in the third and fourth components. The X functions are only defined for 
n a > so components of the n a = wave functions which contain indices outside of 
this range have zero weight on those components. The wave functions are: 



fa 



d a e 



kaV 



\ — e7 -e„ — 9aXn a +l{r a ) \ 



Xn 



V 



Xn a ( r a) 
9aXn a { r a) 



d a e 



( tt>aXn„+l(ra) \ 
UXn a -i{r a ) 

t)aXn a (r a ) 

V laXnS r a) / 



d a 



2(n 



1) 



+ 1 



St 



2n r 



1 



fJa 



(£aS-£ ) 2 

[(£aS + E a ) 2 



with r a = x/X B 



7i £ S + E a 
q a \ B and 

1 



(£ a S + £ a ) 2 
- 2n a l , 



+ 1 



Xn a { r a) 



-.e-^ 2 H na {r a ) 



Irradiated bilayer graphene 



13 



In addition to these functions, there are also two '^-states' with energy ±17/2 which 
account for the doubled degeneracy of the zero energy Landau levels in unbiased 
bilayer graphene. The wave functions of these states are 

/ Xo(r a ) \ 



a ik a y 



<P8a 



V 



Appendix B. Matrix elements of J"(r) 



For convenience, we restate the definition of F(t) = %(t)- 



lK4z = Hn 



-Hu + h( T )- 



ifijfp- The matrix elements contain several terms. The terms relating to the static 
Hamiltonian and the time derivative are 

{{j'X^Ho + Hu-ih^jX)) 



(E x + jhCtA)Sjj'5 x ,x' 



The matrix elements of the term associated with the irradiating field h(r) depend on 
the specific nature of the field and the wave functions. For example, linearly polarized 
light can be described by the vector potential Ai; n (r) = Acos(Qt) [cos 8, sin 8] where 
8 is the angle of polarization in the plane of the graphene with respect to the x axis. 
This yields 



£v F eF 

lllin\T) = — <7 X 







cos^at)- 



n A ~* " V e 

Since we discuss radiation incident at the perpendicular to the graphene plane, this 
angle is immaterial and we can substitute 8 = without loss of generality. On the 
other hand, circularly polarized light is given by A±(t) = A [cos(fl at), ±sin(Q at)] 
where the positive or negative sign corresponds to right- or left-handed orientations 
of the polarization: 

e * iQAT 



£v F eF 
h±{r) = — <j x 







ic field 



Appendix B.l. Zero 

In this case the matrix elements of linearly polarized light (taking 8 = 0) are zero for 
states with unequal wave vectors, and as follows for states in the same valley and with 
the same spin, which also have identical wave vector k: 

((j'ka'\h lin (T)\jka)) = —£=—(5 j >j + i+6 j > J - 1 )C a C a >5 u Jz,t a , x 



e i6k 



fa fa 



■(£5 + E a ,)(Z5 + E a ) 



The matrix elements of states with different spin or valley are zero. Similarly, the 
matrix elements of the Hamilonian associated with circularly polarized radiation are 
diagonal in the wave vector, spin, and valley and are given by 

({j'kX'\h ± (r)\jkX)) = -^LcxC x ,S Ux 8 Uxl x 



X < Sj> 4 



f a f a > (£5 + £ a ,)(£H + £7 a ) 



^ + E a 



7? 
1 

£S + E a , 



+ 



f a f a > (tiZ + E a )(ZZ + E a , 



^-E a 
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Appendix B.2. Quantizing magnetic field 
The matrix elements in this case are 

(0'V|/iiin(T)|ja)) = -j=^ (6j,, j+ i + Sj.j-i) d a d a/ 5 ka , K ,6 u J 6 ,z a , x 



fr , ( 9a9a' 1 



9a9a> 



and 



V2iveF 

n 



(0"a'|/i ± (r)|ja)) = d a d a .6 katka ,6 u J u<t , x 



£5 - E a , £S + £ a 
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